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Abstract 

The recently-discovered narrow charmonium resonance near 3870 MeV is interpreted as a 
hadronic molecule whose constituents are the charm mesons and D*^ or and D*^. Be- 
cause of an accidental fine-tuning of the molecule to very near the D°D*° threshold, it has some 
universal properties that are completely determined by the unnaturally large D^D*^ scattering 
length a. Its narrow width can be explained by the suppression by a factor of 1/a of decay modes 
other than the decay of a constituent D*^ or D*^. Its production rates are also suppressed by a 
factor of 1/a. A particularly predictive mechanism for generating the large scattering length is the 
accidental fine-tuning of a P-wave charmonium state to the D^D*^ threshold. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40. Gx 
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I. INTRODUCTION 



The discovery of the J /ip and other charmonium resonances in 1974 played a crucial role 
in the construction of the gauge theories of the strong and electroweak forces that constitute 
the Standard Model of particle physics. Subsequent studies of the spectroscopy, decays, 
and production of both charmonium and bottomonium resonances have played important 
roles in the development of quantum chromodynamics (QCD), the gauge theory of the 
strong interactions. The recent unexpected discovery by the Belle collaboration of a narrow 
charmonium resonance near 3.87 GeV 0] has presented a new challenge to our understanding 
of QCD. 

The new charmonium state X(3870) was discovered in electron-positron collisions through 
the 5-meson decay -B^ — >■ K^X followed by the decay X — >■ J/ip'K^'K~. Its mass was 
measured to be Mx = 3872.0±0.6±0.5 MeV It is narrow compared to other charmonium 
states above the threshold for decay into two charm mesons: the upper bound on the width is 
Tx < 2.3 MeV. The discovery has been confirmed by the CDF collaboration who observed 
X through J/ipTi^TT~ events inproton-antiproton collisions and measured its mass to be 
Mx = 3871.4 ± 0.7 ± 0.4 MeV 0. 

There have been several recent papers discussing the possible interpretations of the 
X(3870) 0, 0, S S 1^- The most conventional interpretations are previously undis- 
covered states in the charmonium spectrum, such as one of the lowest D-wave states with 
spin/parity /charge-conjugation quantum numbers J^*^ = 2 or 2 ^ or one of the first 
excited P-wave states with J^'^ = 1++ or 1+". A more exotic possibility is a "hybrid char- 
monium" state in which a gluonic mode has been excited. Another possibility, a D^D*^ or 
D^D*^ molecule, is motivated by the fact that the X(3870) is extremely close to the thresh- 
old 3871.2 ±0.7 MeV for decay into the charmed mesons D^D*^. The possibility of hadronic 
molecules formed from charm mesons was suggested long ago 0| . The most favorable 
channels for forming a molecule from the pion-exchange interaction are the P-wave channel 
with J^'-^ = and the S-wave channel with J^*-"' = 1"*"+ j^. 

In this paper, we explore the consequences of identifying X(3870) as an S-wave 
1)0 jj*o ^ jjo £,*o molecule. The tiny binding energy of the molecule implies that the D^D*^ 
scattering length a is unnaturally large. The molecule therefore has properties that depend 
on a but are insensitive to other details of the interactions of and Z)*°, a phenomenon 
called "low-energy universality." In Section IHl we discuss the implications of low-energy 
universality for the wavefunction of X and describe two possible mechanisms for generating 
the large D^D*^ scattering length. In Section IIIIl we discuss the implications of low-energy 
universality for decays of X. One mechanism for generating the large D^D*^ scattering 
length, the fine-tuning of the energy of a P-wave charmonium state to the D^D*^ thresh- 
old, gives a particularly distinctive pattern of branching fractions. Low-energy universality 
gives highly nontrivial predictions for 3-body systems, such as D^D^D*^. Unfortunately, as 
shown in Section HV1 the spectacular possibility of shallow D^D^D*^ molecules called Efimov 
states can be excluded. In Section |3 we present a nonrelativistic effective field theory that 
illustrates the two mechanisms for generating a large D^D*^ scattering length. Our results 
are summarized in Section IVll 



2 



II. LOW-ENERGY UNIVERSALITY 



We will assume that the closeness of Mx to the D^D*^ threshold is no accident and that 
X(3870) is indeed a hadronic molecule whose constituents are the charm mesons and 
Z)*° or and D*^ . What makes this molecule unique among all the hadrons that can be 
interpreted as 2-body bound states of other hadrons is its extremely small binding energy. If 
the low-energy interaction between two hadrons is mediated by pion exchange, the natural 
scale for the binding energy of a molecule composed of the two hadrons is / {2m^ed) ■, where 
mred is their reduced mass. The natural energy scale for a D^D*^ molecule is about 10 MeV. 
The binding energy of the X(3870) (which is positive by definition) has been measured to 
be B2 = —0.5 ± 0.9 MeV. Thus it is likely to be less than 0.4 MeV, which is much smaller 
than the natural energy scale. The only other two-body bound state of hadrons whose 
binding energy is known to be small compared to the natural energy scale is the deuteron. 
Its binding energy is 2.4 MeV, which is small compared to the natural energy scale of 20 
MeV for a pn molecule. 

We will further assume that X(3870) is an S-wave bound state of D^D*^ or D^D*^, be- 
cause this has particularly interesting implications. Since the constituents have quantum 
numbers 0~ and 1~, the J^*" quantum numbers of the molecule must be 1"^^ or 1+^. The 
interaction between and Z)*" at energies less than m^/mD ~ 10 MeV is dominated by 
the S-wave channel and can be described by a single parameter: the S-wave D^D*^ scat- 
tering length a. A shallow S-wave bound state implies an S-wave scattering length that 
is large compared to the natural length scale 1/m^. The low-energy few-body observables 
for nonrelativistic particles with short-range interactions and a large scattering length have 
universal features that are insensitive to the details of the mechanism that generates the 
large scattering length. This phenomenon is called low- energy universality. If a > 0, the 
simplest universal prediction is that there is a shallow 2-body bound state whose binding 
energy B2 for sufficiently large a approaches 

^ 2mrcdO^ ' ^ ^ 

where m^ed is the reduced mass of the two constituents. If the binding energy of the X(3870) 
were measured, the D^D*^ scattering length a could be predicted using (jT)) with mj-ed = 
m£)om£)*o /{rrLfjo + m/5.0). For example, if the binding energy of X(3870) were 0.5 MeV or 
0.1 MeV, the scattering length would be 6.3 fm or 14.2 fm, respectively. These are both 
much larger than the natural length scale 1/^1^^ = 1.5 fm. 

Low-energy universality has other implications for the interpretation of X(3870) as a 
D^D*^/D^D*^ molecule. There is a universal prediction for the D^D*^ or D^D*^ wavefunc- 
tion: 



(2) 



Voloshin has exploited this universal wavefunction to calculate the momentum distributions 
for the decays X — > D^D^tc^ and X D^D^'-f [6J. There are also components of the 
wavefunction that correspond to other hadronic states with the same J^*" quantum numbers. 
If J^'^ = 1+"*", they include the P-wave charmonium states Xci(l-P) and Xci(2-P). If J^'^ = 
1+^, they include the P-wave charmonium states hc{lP) and hc{2P). In either case, they 
also include D^D*^ /D^D*^ states and D*^D*^ states. In an appropriate hadronic basis, 
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the quantum state for X(3870) can be written 

+ J2ZT\H). (3) 



H 



where ip{p) is the Fourier transform of the D^D*^/D^D*^ wavefunction and the sign ± 
is determined by the charge conjugation quantum number (7 = ±1. The other hadronic 
states H can be discrete states, such as Xci(2-P) or hc{2P), or continuum states, such as 
D^{p)D*~{—p). The probabihty factors Zj^d* and are real and positive, and they 
add up to 1. Low-energy universahty imphes that as the scattering length a increases, the 
probabilities for states other than D°D*° or D^D*^ decrease as 1/a and the D°Z)*°/Z)°L)*° 
wavefunction approaches (j21). In the limit a —>■ oo, the state becomes a pure D^D*^/D^D*^ 
molecule. 

A scattering length that is large compared to the natural length scale necessarily requires 
a fine-tuning. In the case of the D^D*^ molecule, the fine-tuning parameters can be iden- 
tified with the up and down quark masses m„ and m^. The masses of and D*^ are 
sensitive to m„, because these hadrons contain an up quark as a constituent. The D^D*^ 
potential is sensitive to m„ and irid through the pion mass. There are two distinct mecha- 
nisms for generating a large D^D*^ scattering length. The first mechanism is a fine-tuning 
of parameters that have a large effect on the D^D*^/D^D*^ channel without significantly 
affecting other channels. This could be a fine-tuning of the range and depth of the D^D*^ 
potential so that there is a bound state very close to threshold and thus a large scattering 
length. Equivalently, it could be a fine-tuning of the masses of the and D*^ to obtain 
a bound state very close to threshold in the D^D*^ potential. This mechanism requires 
the quantum numbers of X to be J^^ = because this is the only S-wave channel for 
which the potential due to pion exchange is sufficiently attractive to produce a bound state 
[3]. In the limit a — > oo, the probabilities for components of the wavefunction other than 
£)0£)*o 1)0 D*o g(.g^ig rpj^jg illustrated in Section |V| using an explicit field 

theory model. If the energy gap uh between the state H in ^ and the D^D*^ threshold 
is much greater than the natural energy scale m^/mred, a more complete estimate of the 
dimensionless suppression factor in Zh is / (m'^^^i'jja) . For most channels, the energy gap 
uh is much larger than the natural energy scale. For example, the energy gap for Xcii^P) is 
uh = —360 MeV. If am^r ^ 1, the suppression factor ml / (m'^^^i/jfa) is 1/1600 for B2 = 0.5 
MeV and 1/3500 for B2 = 0.1 MeV. If this is the correct mechanism for generating the large 
scattering length, we can probably neglect all components of the wavefunction other than 
^0^*0/^0^*0 ^-^^ gg^ 2dd* ~ 1. 

A second mechanism for a large D^D*^ scattering length is an accidental fine-tuning of one 
of the P-wave charmonium states Xci(2-P) or hc{2P) to the D^D*^ threshold. X(3870) will 
have the same J^^ quantum numbers as the charmonium state: 1"'""'" in the case of Xci(2P) 
and l"*"" in the case of hc{2P). This mechanism is analogous to the Feshbach resonances [11| 
that can be used to control the scattering lengths for atoms by adjusting the magnetic field 
jl^ [13L l3 |. Feshbach resonances are currently being used to tune the scattering lengths 



for atoms to arbitrarily large values in order to study Bose-Einstein condensates of bosonic 
atoms and degenerate gases of fermionic atoms in the strongly-interacting regime. In the 
case of the D^D*^ molecule, the fine-tuning parameter can be identified as rriu, which can 
shift the and Z)*° masses, thus changing the energy gap z/ between the Xci{2P) or hc{2P) 
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and the D^D*^ threshold. In potential models, which ignore the coupling of charmonium 
states to continuum channels such as DD and DD*, the estimates of the energy gap for 
Xci(2P) or hc{2P) are both u ^ 90 MeV jisl . EgI ]. The predictions of these models for the 
mass of the ID state ■?/'(3770) with J^*" = 1 are too large by about 50 MeV, so they 
are also likely to overpredict the masses of the Xci(2-P) and hc{2P). The error presumably 
arises mostly from the neglect of coupled-channel effects, which are sensitive to the light 
quark masses. If the coupled-channel effects shift the 2P charmonium states down by about 
90 MeV relative to the DD* threshold, they could fortuitously tune the energy gap u for 
Xci(2P) or hc{2P) to be smaller than the natural low-energy scale m^/mi:) ^ 10 MeV 
associated with pion exchange. In this case, a resonant interaction between the Xci(2P) or 
hc{2P) and D^D*^ states generates a large D^D*^ scattering length a that increases as l/u. 
If a > 0, there is a shallow bound state whose binding energy approaches ([T]) as — *• 0. 
In the expression Q for the quantum state, the hadrons should be interpreted as those in 
the absence of the fine-tuning that generates the resonant interaction. In the limit a ^ oo, 
the probability for Xci{2P) or for hc{2P) scales as 1/a, as do the probabilities Zh 
for all other channels besides D^D*^/D^D*^. This will be illustrated in Section IVI using an 
explicit field theory model. The probability Z^ or Zh includes a dimensionless suppression 
factor l/(am7r), whose value is about 1/4.3 for B2 = 0.5 MeV and 1/9.7 for B2 = 0.1 MeV. 
There may also be further suppression from a small numerical coefficient associated with 
Zweig's rule, because the processes Xci(2P) D^D*^ and hc{2P) — >• D^D*^ require the 
creation of a light quark-antiquark pair. If this is the correct mechanism for generating the 
large scattering length, we can probably neglect all components of the wavefunction other 
than D^D*^/D^D*^ and Xci(2P) or K{2P). We can then set Zdd* ^ I - Z^in the case 
jpc ^ i++ g^^^ ^ ^ 1 - Zft in the case J^^ = l+~. 



III. DECAYS 

An important requirement for any interpretation of the X(3870) is that it provide an 
explanation for its narrow width. The upper bound Tx < 2.3 MeV implies that the width 
of X is more than an order of magnitude smaller than that of the D-wave state ip{3770). 
According to our interpretation, X is below the D^D*^ threshold and it therefore cannot 
decay into D^D*^. Its quantum numbers = 1^ forbid a decay into D^D^ or D^D~. 
It can however decay into D^D^ir^ or D^D^'-f by the decay of a constituent D*^ or Z)*° of 
the molecule. It can also decay into a lighter charmonium state by a radiative or hadronic 
transition or into light hadrons via a process in which the charm quark and antiquark 
annihilate. 

We first consider the contribution to the width Tx from the decay of a constituent 
D*^ into D^n^ or D^^. The width r[D*^ —>■ D^n^] can be deduced from the measured 
width of the D*~^, the branching fraction for D*~^ D^n^ or D*^ — D^tt^, and isospin 
symmetry: r[D*^ D'^n^] = 31 ± 7 keV. We have treated the difference between the 
branching fraction for D*^ — > D^tt^ and twice the branching fraction for D*~^ D^tt^ as a 
systematic error associated with isospin symmetry breaking. The width T[D*^ — * D^^] can 
then be deduced from the measured branching fractions for D*^ D^n^ and D*^ D^'-j: 
T[D*^ —* D^'-f] = 19 ± 5 keV. The contributions to the width Tx from these decays are 

r[X Z}°D°7r°] = War[D*° ^ D\% (4a) 
r[X^D°D°7] = Zdd*C^T[D*'' ^ D^'-f]. (4b) 
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The factors Ct^ and take into account interference between the decay of D*^ from the 
]jO£)*o component of the wavefunction and the decay of D*^ from the D^D*^ component. If 
the charge conjugation quantum number of X is C = +1, there is constructive interference in 
the decay X D^D^n^ and destructive interference in the decay X D^D^'-). If C = — 1, 
the pattern is reversed 0]. If C = +1 (or —1), the coefficient C,^ ranges from about 1.5 (or 
0.13) if the binding energy B2 is 0.5 MeV to about 2.2 (or 0.56) if B2 = 0.1 MeV and to 2 
if B2 = 0. The coefficient ranges from about 0.58 (or 3.42) if B2 = 0.5 MeV to about 
1.36 (or 2.64) if B2 = 0.1 MeV and to 2 if ^2 = 0- The probability factor Zdd* is close 
to 1. Thus the lower bound on the width provided by the sum of Paj) and ()4b|) ranges from 
56 ± 11 keV (or 69 ± 16 keV) if B2 = 0.5 MeV to 92 ± 17 keV (or 67 ± 13 keV) if ^2 = 0.1 
MeV to 99 ± 17 keV if B2 = 0. 

There can also be decays that proceed through the decay of a constituent D*^ or D*~ 
from the D~D*^ or D^D*^ component of the wavefunction. The decays of X into D^D^n^, 
D^D'TT^, or D^D^n^ are forbidden by phase space. The decay X — > D~^D~'~f is allowed, 
but it is suppressed by the small decay rate r[D*^ D^l] = 1-5 ± 0.5 keV and also by 
the small probability for the D~D*^ /D^D*~ component of the wavefunction. Thus the 
contribution to the width Tx from the D~D*^ / D^D*~ component of the wavefunction can 
be neglected. 

The remaining decay channels of the X(3870) are radiative transitions to lower charmo- 
nium states such as X — > ip{2S)'y or X — > ?7c(2S')7, hadronic transitions to lower charmonium 
states such as the discovery mode X — > J/ip7r~^n~ , and annihilation decays into light hadrons 
such as X — ^ pp. These decays proceed through the short- distance part of the D^D*^ / D^D*^ 
wavefunction or through other components of the wavefunction, such as Xci(2-P) or hc(2P). 
All these contributions to the decay rate are suppressed by a factor of 1/a and go to as 
the binding energy of X goes to 0. 

If the large D^D*^ scattering length arises from an accidental fine-tuning within the 
D^D*^/D^D*^ channel with J^'-^ = 1++, the probability Zh for a component of the wave- 
function with a large energy gap {uhI ^ m^/m-^eA is suppressed by / {mf^^ufja) . Thus 
the radiative transitions, hadronic transitions, and annihilation decays of X(3870) are dom- 
inated by the D^D*^/D^D*^ component of the wavefunction. For example, the amplitudes 
for radiative and hadronic transitions to J/ ip can be approximated by 

^[X-.J/^ + 7] ^ Z'JlJ ^^(p)v^^[D«(p)5*°(-p)^ J/^ + 7], (5a) 

A[X ^ J/^ + h]^ Z]il I ^J(p)V2A[D%p)D*'{-p) J/^ + h], (5b) 

where h consists of light hadrons. An example of a light hadronic state h is the discovery 
channel tt+tt". The light hadronic state h must be odd under charge conjugation. Since D° 
and D*^ have isospin ^ and J/ip has isospin 0, h can have isospin or 1. For p <^ m^r, the 
momentum space wavefunction ipip) can be approximated by its universality limit: 

The momentum integrals in (0) are cut off at large momentum by the p-dependence of the 
transition amplitude for D^D*^ — >■ J/ip + h. For this transition to occur, the heavy c and c 
in the J/ijj must be in the same momentum states as the c and c in the and Z)*°. Thus 
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the transition amplitude includes a factor of the momentum space wavefunction i/jj/^pip) for 
the J/ip. The momentum scale associated with this wavefunction is rricV ~ 700 MeV, where 
V is the typical velocity of the charm quark in chaxmoriium. The amplitude on the right side 
of (0) therefore includes the overlap factor J d^p^{p) ipj/^ipip). The explicit factor of a~^/^ in 
the universal wavefunction in Q combines with a factor of (mcf)"^/^ from the integral to 
give a dimensionless suppression factor. Thus the decay rate scales as l/{amcv) as a ^ oo. 
If a phenomenological framework for calculating radiative and hadronic transitions of D^D*^ 
to quarkonium were available, the rates for radiative and hadronic transitions of X could be 
calculated using equations analogous to (0). 

If the large D^D*^ scattering length arises from an accidental fine-tuning of the P-wave 
charmonium state Xci(2-P) to the D^D*^ threshold, the radiative and hadronic transitions 
and the annihilation decays of X(3870) can also proceed through the Xci(2P) component of 
the wavefunction. These contributions to the decay rate are suppressed by the probability 
factor Zy., which scales as l/{am.,^). Although this suppression factor has the same power of 
a as in the D^D*^ contribution, the Xci(2-P) contribution may be numerically larger because 
of a factor of mcv/m.,^ in the ratio of the suppression factors. Thus decays of X into modes 
that are possible final states of the decay of Xci(2-P) are likely to be dominated by the 
Xci(2P) component of the wavefunction. The rate for these decays will be given by the rate 
for the corresponding decays of the Xci(2P) in the absence of the fine-tuning multiplied by 
the probability factor Z-^. For example, the decay rates for radiative transitions to J /ip and 
for hadronic transitions to J/ip via the emission of a light hadronic state with total isospin 
are 

r[X ^ J/V^ + 7] ^ Z^V[xci{2P) ^ J/V^ + 7], (7a) 
V[X^ J/^ + /i,=o] ~ Z^V[xci{2P) ^ J/^ + /i,=o]. (rb) 

The last factors in (fTaj) and (ffE]) are the decay rates for Xci(2-P) assuming that it has mass 
mx but ignoring the resonant interaction with D^D*^ . Hadronic transitions in which the 
light hadronic state h has total isospin 1, such a.s X ^ J /ip + cannot proceed through 
the Xci(2P) component of the wavefunction. They must therefore be dominated by the 
D^D*^/D^D*^ component. The amplitude for such a transition can be approximated by 
fl5b|) . There are well-develop ed p henomenological frameworks for calculating the radiative 
and hadronic transition rates 0,^^ for charmonium states such as Xci(2-P). If the rates for 
several radiative transitions or hadronic transitions with isospin were measured and found 
to be all suppressed relative to the predictions for XciC^P) decays by a common factor Z^, it 
would be strong evidence in favor of this fine-tuning mechanism. The total width of Xcii'^P) 
in the absence of the resonant interaction with D^D*^ must be significantly larger than the 
width of Xci(l-P), which is about 1 MeV, since Xci(2P) has more decay channels and the 
decays have larger phase space. However the suppression from the probability factor Z^ 
could reduce this contribution to the width Fx so that it is comparable to or even smaller 
than the contribution from the decay of D*^ or Z)*°. 

If the large D^D*^ scattering length arises from an accidental fine-tuning of the P-wave 
charmonium state hd^P) to the D^D*^ threshold, the radiative and hadronic transitions 
and the annihilation decays of X can also proceed through the hc{2P) component of the 
wavefunction. The decay rates for radiative transitions and for hadronic transitions with 
total isospin would be given by the corresponding decay rates of hc{2P) multiplied by a 
probability factor Zh- 

The identification of X(3870) as a shallow S-wave molecule also has implications for its 
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production rate in high-energy colhsions. As the D^D*'^ scattering length a increases, the 
production rate decreases as 1/a. If the large value of a arises from an accidental fine-tuning 
within the D^D*^ channel, the production will proceed primarily through the creation of 
and Z)*° (or and with small relative momentum of order 1/a. In this case, the 
suppression factor of 1/a in the production rate comes from the factor a~^/^ in the universal 
wavefunction If the large value of a arises from an accidental fine-tuning of the P-wave 
charmonium state Xci(2-P) to the D^D*^ threshold, the production rate may be dominated 
by production of the Xci(2P). In this case, the suppression factor of 1/a comes from the 
probability for the Xci(2P) component of X. Similarly, if the large value of a arises 
from an accidental fine-tuning of hc{2P) to the D^D*^ threshold, the production rate is 
suppressed by a factor of 1/a from the probability Zh for the hc{2P) component of X. 



IV. ABSENCE OF EFIMOV STATES 

The most remarkable predictions of low-energy universality, which were discovered by 
Efimov [3], occur in the 3-body sector. At sufficiently low energies, the effective inter- 
action between three nonrelativistic particles with short-range forces can be described by 
an effective potential Kff(-R) that depends only on the hyperspherical radius R, which is a 
weighted average of the separations of the three particles |2fl]. If the scattering length is 
large compared to the range i of the force, the effective potential in the region i <^ R <ti \a\ 
is scale-invariant. In the case of identical particles of mass m, the hyperspherical radius is 
just the root-mean-square separation of the three pairs of particles and the scale-invariant 
potential is 

where Aq is the minimum of the nontrivial solutions to 

\/3A^/2 cos(7rA^/V2) = 8 sin(7rA^/V6). (9) 

The minimum solution is Aq = — Sg, where sq ~ 1.00624. In the resonant limit a ^ oo 
in which the scattering length is tuned to be infinitely large, the 2-body bound state has 
zero binding energy and there are infinitely many arbitrarily-shallow 3-body bound states 
called Efimov states. If the particles are identical bosons, the ratio of the binding energies 
of adjacent states approaches a universal constant e'^'^^^° ~ 515.03. The 3-body spectrum in 
the resonant limit has an asymptotic discrete scaling symmetry with discrete scaling factor 
g7r/so ~ 22.7. This symmetry is related to an infrared renormalization group limit cycle [19\ . 
A limit that is more relevant to a physical problem with a large but finite scattering length 
is the scaling limit defined by A — oo with a fixed, where A is the natural momentum scale 
set by the inverse of the range of the interaction. In the scaling limit, the binding energies 
i?3 and -Bg of the shallowest and next-to-shallowest Efimov states for identical bosons are 
in the intervals B2 < < 6.75^2 and 6.75^2 < B'^ < 1406^2, where B2 = l/{ma'^) is the 
2-body binding energy |21,]. If these binding energies are smaller than the natural energy 
scale A^/m, these Efimov states should appear as real states in the spectrum. Thus, there 
should be at least one Efimov state if A^/m > 6.75B2 and at least two if A^/m > 1406i?2- 



As an illustration, we consider helium atoms, which have a large scattering length [2l[ . The 



helium dimer is very shallow: its binding energy B2 ~ 1.3 mK is smaller than the natural 
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low-energy scale h?' /m ^ 400 mK by about a factor of 300. Thus we would expect either 
one or two Efimov states. There are in fact two helium trimers: a ground state and an 
excited state with binding energies B3 ~ 130 mK and B3 ^ 2 mK. Both can be interpreted 
as Efimov states jiH ]. 

The large D^D*^ scattering length raises the exciting possibility of shallow D^D^D*^ 
molecules within 10 MeV of the D^D^D*^ threshold generated by the Efimov effect. Unfor- 
tunately, this possibility can be excluded. The D^D^D*^ sector involves only two identical 
bosons and only two of the three pairs of particles have a resonant interaction with a large 
scattering length. Furthermore a D^D*^ pair can fluctuate into a D^D*^ pair, and the other 

has no resonant interaction with this component of the wavefunction. Low-energy inter- 
actions in the 3-body sector can again be described by an effective potential which in the 
region -C -R <^ \a\ has the scale-invariant form (jHJ. The form of the potential can be 
derived from results given in Ref. ji^]. If we ignore the 8% mass difference between the 
and -D*", the only difference is that the equation for Aq is 

VSA^/' cos(7rA^/V2) = 2 sin(7rA^/V6). (10) 

Of the factor of 4 difference with Q, one factor of 2 comes from there being only two 
identical bosons instead of three and the other factor of 2 comes from the 3-body system 
being a superposition of a D^D^D*^ molecule and a D^D^D*^ molecule. The minimum 
nontrivial solution to (fTUjl is Aq ~ 0.3533. Since this is positive, the Efimov effect does not 
arise and there are no shallow 3-body bound states. 



V. A FIELD THEORY MODEL 

If we consider only momenta small compared to the natural momentum scale m^r, hadrons 
such as and D*^ can be treated as point particles with pointlike interactions and can 
therefore be described by a local nonrelativistic field theory. In Section |Tll we discussed two 
fine-tuning mechanisms that can generate a large scattering length for D^D*^ /D^D*^ . Only 
one of these mechanism is capable of producing a large scattering length in the 1^ channel, 
but either one is capable of producing a large scattering length in the l^"*" channel. We 
will therefore focus on the possibility J^'^ = for the quantum numbers of X(3870). A 
model that can describe either of the two fine-tuning mechanisms is a nonrelativistic field 
theory with local fields D, Z), D , D, and x for the D*^, D*^, and Xci(2P). The 

hamiltonian density is the sum of mass terms, kinetic terms, and interaction terms: 

^mass = rriDO {D^ D + D^D) + mz5*o (Dt ■ D + • D) 

+ {17100 + mD*o + uo)x^ ■ X (11a) 
Hkin = -^m^l {D^V^D + D^V^D) - \m-J,, (D^ • V'D + ■ V'D) 

-i(mBO + mB,o)-ix^- V^x, (lib) 
Hint = Ao {Dti + DD)^ • (DD + DH) 

+go [x'-{Dti + DD) + (DD + DTi)^-x] , (He) 

where Aq, Qq-, and z/q are bare parameters that require renormalization. A similar field theory 
has been used to describe the behavior of cold atoms near a Feshbach resonance i^l . If we 



impose an ultraviolet cutoff A on loop momenta and drop terms that decrease as inverse 
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powers of A, the cutoff dependence can be removed by eliminating Aq, go, and z/q in favor of 
renormalized parameters X, g, and u defined by 



A = Z-'\o, (12a) 
g = Z^'go, (12b) 
u = iyo + [Z~' - l]gl/Xo, (12c) 



2"^rcdA. (14) 



where the renormahzation constant Zx is 

2 

Zx = l + ^mredAoA (13) 

and mj-ed = Tnuomij*o/[m£,o + mD*o) is the reduced mass. Note that the combinations 
go/Xo = g/X and i^o ~ go/ — ^ 9"^/^ are renormahzation invariants. 

The natural scale for the ultraviolet momentum cutoff is A ~ m^r. The natural 
magnitude for the bare coupling constant Aq can be deduced by dimensional analysis: 
|Ao| ~ 1 / {mj-cdrriTr) . This can be made evident by writing the renormahzation condition 
(|12a|) in the form 

1 _ 1 2 
A Ao vr 

If the renormalized coupling constant A is fixed and A is sufficiently large, Aq must scale 
like (mi-edA)^^ to compensate for the effect of the ultraviolet cutoff. The natural magnitude 

for go is (^ml/'^ /rrij-cd, where the factor of /m^^d comes from dimensional analysis and 
C is a numerical suppression factor associated with the violation of Zweig's rule by the 
process Xci(2-P) — > D^D*^, which requires the creation of a light quark- ant iquark pair. The 
renormahzation condition ()12b|) implies that the numerical suppression factor ( is stable 
under renormahzation and does not require fine-tuning. There is no natural magnitude 
for the bare parameter z/q: it is completely adjustable. In the absence of fine-tuning, the 
renormahzation constant Zx in ()13|) is comparable to 1. The renormahzation conditions 
fll2ap . ()12b|) . and ()12c|) then imply that the natural magnitudes of the renormalized coupling 

constants are |A| ~ (mred'^Tr)"^, \g\ ~ C'lTh/'^ /mj-ed, and |z/| ~ max(|z/o|, C^^^^/^^red)- 

The 2-body scattering amplitude in this model can be calculated analytically. The scat- 
tering length is 

"^red A g^ 



TV \ V 



• (15) 



The natural magnitude for |a| is XjniT^. The scattering length can be made unnaturally large 
either by making A sufficiently large, which corresponds to tuning the potential between 
D° and Z)*°, or by making v sufficiently small, which corresponds to tuning the energy 
gap between the Xci(2-P) and the D^D*^ threshold. In either case, low-energy universality 
implies that as a increases, the binding energy of the molecule approaches and the D^D*^ 
or D^D*^ wavefunction approaches (j21). 

The first mechanism for generating a large scattering length is to make A unnaturally 
large: |A| ^ |Ao|. This can be accomplished by tuning Aq towards the critical value 
— 7r^/(2mi.edA), so that there is a near cancellation between the two terms on the right 
side of ()14|). This fine-tuning makes the renormahzation constant Zx much less than 1. The 
renormahzation condition ()12b|) implies that this fine-tuning also increases the strength of 
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the effective coupling constant between x ^i-nd DD: \g\ ^ l^fol- This is also evident from the 
fact that g/X = qo/Xq is a renormalization invariant. There is a limit to how large the scat- 
tering length can be made using this mechanism. When Z\ becomes smaller than (y'o/l-^o^'ol, 
the Qq/Xq term in (|12c|) begins to dominate over the vq term. In this case, both terms in 
the scattering length (|15p become large and they tend to cancel each other. Thus, with this 
mechanism, the maximum magnitude of the scattering length is of order (Ao/5'o)^'^rcd|t'o| 
which is of order C~^'^red|t'o|/''^^- 

The second mechanism for generating a large scattering length is to make v sufficiently 
small. This can be accomplished by tuning uq towards the critical value —[Z^^ — IJ^^q/Ao 
for which there is a near cancellation between the two terms on the right side of p2c|) . The 
scattering length can be made arbitrarily large using this mechanism. 

The calculation of the binding energy B2 of X and of the probability Z-^ for the Xci (2-P) 
component of the wavefunction can both be reduced to the solution of a cubic polynomial. 
The binding momentum k defined by B2 = k."^/ {2m^cd) satisfies the cubic equation 

+ 2mj.cdi^ = — —Xk [k^ + 2mi.ed — ■ (16) 

In either of the two limits A ^ 00 or ^ 0, one of the three roots of this polynomial has 
the limiting behavior k, 1/a. If a > 0, the probability Z^ for the Xci{2P) component of 
the wavefunction is 

^ 2n\[K^ + 2m,^d{jy-9yX)]^ ATtKj + 2m,^4u - g^X)' ^ ^ 

After expressing the observables B2 and Z^^ as functions of a and the renormalization in- 
variants g/X and u — g'^/X, they can be expanded in powers of 1/a: 

B2 ^ 1 (i <9l\f ^ \^ ^^g^) 

^ 2mreda?\ m^^^{u - g'^/X^a ' ' / ' 

. (18b) 

For any fine-tuning that produces a large scattering length, the bare coupling constants 
approach definite limiting values and therefore the renormalization invariants g/X and u — 
g^/X approach definite limiting values. Thus the probability decreases like 1/a. This 
illustrates our assertion that with the probability for states other than D*^D*^ or f)*^D*^ 
scales as 1/a. 

We proceed to discuss how the decays of X(3870) could be described within this effective 
field theory. In the decay X — >■ D^D^ti^, which is dominated by the decay of a constituent 
£,*o ^*o^ ^Yie typical momentum of the final or is 40 MeV, which is much smaller 
than the natural momentum scale m^r. Thus this decay can be described within the effective 
theory by introducing a vr*^ field into the lagrangian with an interaction term that allows the 
decay Z)*" — D^n^. In the decay X — >■ D^D^'y, which is also dominated by the decay of a 
constituent D*^ or D*^, the typical momentum of the recoiling or Z)° is 140 MeV, which 
is comparable to the natural momentum scale m^r. Thus this decay need not be described 
accurately within an effective theory in which hadrons are treated as point particles with 
pointlike interactions. 

The radiative and hadronic transitions and the annihilation decays of X(3870) produce 
particles with momenta larger than the rriTr- They therefore cannot be described explicitly 
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within an effective theory in which hadrons are treated as point particles with pointhlce 
interactions. The inclusive effects of these decays can however be taken into account im- 
plicitly through local terms in the hamiltonian density. The inclusive effects of transitions 
of D^D*^ or D^D*^ to charmonium states and of their annihilation into light hadrons can 
be taken into account through an imaginary part of the bare coupling constant Aq. The 
inclusive effects of transitions of Xci (2-P) to other charmonium states and of its decays into 
light hadrons can be taken into account through an imaginary part of the bare parameter 
z/q: Imz/Q = — |r^^-^(2P). The imaginary part of can take into account interference effects 
associated with transitions of D^D*^ and Xcii'^P) to the same final states. If the parameters 
Ao, go, and z/q have small imaginary parts, the scattering length (fT3j) is complex-valued with 
a small imaginary part. If a fine-tuning makes the real part of a large, the binding energy of 
X is given by the real part of the expression (0). Its imaginary part should be interpreted 
as ^ATx, where ATx is the contribution to the width from the decays whose effects are 
taken into account through ImAo, Im(yfo, and Imuo- At first order in the imaginary parts of 
Ao, go, and z/q, the contribution to the width is 

(19) 



ATx = — 



2aA/7r)2(-ImAo) + 2-(l - 2aA/7r)Im^o + ^(-Imz/Q 



If we express gjv iw terms of a and the renormalization invariants, we see that it increases 
linearly with a: gjv = a{g/X)/{v — g'^/X)- Thus all three terms in (jl9p scale as 1/a in the 
limit a ^ oo. This scaling behavior is in agreement with that deduced in Section UTTl 



VI. SUMMARY 

We have explored the implications of low-energy universality for the identification of 
X(3870) as an S-wave D^D*^/D^D*^ molecule. Its shallow binding energy requires some 
fine-tuning mechanism to generate a large D^D*^ scattering length a. Two possible mecha- 
nisms are an accidental fine-tuning of parameters associated with the D^D*^ sector and an 
accidental fine-tuning of a P-wave charmonium state to the D^D*^ threshold. A field theory 
model that illustrates both of these mechanisms was presented. With either mechanism, the 
probabilities for components of the wavefunction other than D^D*^ or D^D*^ are suppressed 
by a factor of 1/a. The decay rates into modes other than those associated with decay of a 
constituent D*^ or D*^ are also suppressed by a factor of 1/a. The assumption that the large 
scattering length arises from the fine-tuning of Xci(2P) or hc{2P) to the D^D*^ threshold 
is particularly predictive. The decay rates for radiative transitions and for hadronic transi- 
tions via emission of light hadrons with total isospin should differ from the corresponding 
decay rates of Xci(2-P) or hc{2P) in the absence of the fine-tuning by a common suppression 
factor. Low-energy universality also has nontrivial predictions for 3-body systems, such as 
D'^D^D*^, although the spectacular possibility of Efimov states can be excluded. In conclu- 
sion, if the X(3870) is indeed an S-wave D^D*^/D^D*^ molecule, it will provide a beautiful 
example of the remarkable phenomenon of low-energy universality. 

This research was supported in part by the Department of Energy under grant DE-FG02- 
91-ER4069. 
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